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Open and Hidden Charm Production in 920 GeV Proton-Nucleus Collisions 

U. Husemann 2 for the HERA-B Collaboration 
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The HERA-B collaboration has studied the production of charmonium and open charm states in collisions of 920 GeV 
protons with wire targets of different materials. The acceptance of the HERA-B spectrometer covers negative values of xp up 
to xp = —0.3 and a broad range in transverse momentum from 0.0 to 4.8GeV/c. The studies presented in this paper include 
J/y differential distributions and the suppression of J/i// production in nuclear media. Furthermore, production cross sections 
and cross section ratios for open charm mesons are discussed. 



1. Introduction 

Studying charmonium and open charm hadropro- 
duction provides a good test of the theoretical produc- 
tion models available and their modifications in nu- 
clear media. Since many of these models rely on ex- 
perimental data to adjust their free parameters, a com- 
prehensive test of the model predictions in a broad 
kinematic range and with good precision is desirable. 

HERA-B is in the position to extend previous mea- 
surements of J/V production and nuclear suppression 
to negative xp values. D° and D + production cross 
sections and cross section ratios are determined with 
high accuracy. 

2. The HERA-B Detector and Trigger 

The HERA-B detector is a fixed-target spectrom- 
eter with large angular acceptance. The detector in- 
cludes subdetectors for vertexing, tracking, and parti- 
cle identification, as shown in Fig.[0 

Protons with energies of 920 GeV are brought into 
collisions with an internal wire target. The target con- 
sists of two stations of four wires each, separated by 
4cm along the beam axis. Each wire can be moved in- 
dependently in the halo of the proton beam to adjust 
the interaction rate. 

The vertex detector system (VDS) allows for re- 
construction and separation of primary and secondary 
vertices. It consists of eight double-layers of double- 
sided silicon strip detectors. The tracking detector is 
divided into an inner part consisting of micro-strip 
gaseous chambers with gas electron multiplier foils 
and an outer part of honeycomb drift chambers. 

The particle identification devices include a ring- 
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Figure 1 . Side view of the HERA-B detector. 



imaging Cerenkov counter (RICH), an electromag- 
netic calorimeter (ECAL), and a four-layer muon de- 
tector. 

In HERA-B, a multilevel trigger system is em- 
ployed to enrich lepton pairs from i/y decays. Pre- 
triggers in the muon detector and ECAL provide track 
seeds for the first level trigger (FLT). The FLT per- 
forms a track search in the outer tracker. In the sec- 
ond trigger level, the tracking of the FLT is confirmed, 
tracks are extrapolated to the VDS, and a vertex fit is 
performed for pairs of leptons. Accepted events are 
reconstructed online. Alternatively, a minimum bias 
trigger can be utilized, requiring either a minimum 
number of photons in the RICH or a minimum energy 
deposit in the ECAL. 

3. The 2002/2003 Data-Taking Period 

The results presented in this paper are based on 
data taken with the HERA-B detector from Novem- 
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Figure 2. J/y signal in J/y — > decays (left) and 

from J/y — > e + e~ decays (right), extracted from min- 
imum bias triggered data. 

ber 2002 to February 2003. In dilepton triggered 
data, about 170,000 J/y — ► p + p~ decays and about 
150,000 J/y - ¥ e + e~ decays have been reconstructed. 
In addition, about 200 million events have been 
recorded with the minimum bias trigger. 

4. Charmonium Production 

4.1. J/y Cross Section 

A clean J/y signal is observed in the data recorded 
with the minimum bias trigger. The invariant mass 
distributions for the dimuon and the dielectron chan- 
nels are shown in Fig. [2] From these signals the J/y 
production cross section can be extracted without trig- 
ger bias. 

4.2. J/y Differential Distributions 

From data recorded with the dilepton trigger, differ- 
ential cross sections of J/y production are determined. 
As an example, in Fig. [5] the J/y — > e + e~ yield for 
25% of the full statistics is shown as a function of 
the transverse momentum pj. The HERA-B detector 
covers a broad kinematic range of 0.0 to 4.8 GeV/c. 
The pj spectrum is fitted with the parametrization 
function 
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in which the normalization constant A and the aver- 
age transverse momentum (pj) are free parameters. 
In Tabled preliminary HERA-B results on (pj) are 
shown. The comparison of different target materi- 
als and with results of the experiments E771 [ 1 1 and 
E789 |2| shows an increase of (pj) with the atomic 
mass number. 
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Figure 3. Transverse momentum distribution 
AN /dp t 2 of J/y — * e + e~ decays. 



4.3. A-Dependence of J/y Production 

Nuclear effects in heavy quark production are com- 
monly parametrized by the power law 



= a pN -A a( - Xe '^ 



(2) 



Here, g p a is the proton-nucleus cross section, g p n is 
the proton-nucleon cross section, and A is the atomic 
mass number of the target. A value of (x(xf,pt) < 1 
corresponds to suppression of heavy quark produc- 
tion. The sources of nuclear suppression in J/y pro- 
duction are either interactions of the partons in the ini- 
tial state with the nuclear medium, or the absorption 
of the cc pre-state or the final J/y, see e.g. Q. 

In HERA-B, the suppression parameter a is ex- 
tracted from data samples taken with a carbon and a 
tungsten target simultaneously. This setup is utilized 
to minimize systematic uncertainties in the detector 
and trigger performance. Using a = N / (Jzfe), Eq. (|2jl 
can be solved for a: 
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The A-dependence measurement hence includes 
measuring three ratios, the ratio of J/y yields, 
Nw/Nq, the ratio of luminosities per wire, ^fc/^w, 
and the ratio of efficiencies, ec/ew- The luminosity 
ratio is extracted from the number of interactions in 
minimum bias events recorded in parallel to the trig- 
gered data. The ratio of efficiencies is determined 
from a detailed Monte Carlo (MC) simulation of the 
HERA-B detector and trigger. Preliminary results on 
a(xp), based on 10% of the full statistics, are pre- 
sented in Fig. |4] Data from different combinations 
of carbon and tungsten wires show results consistent 
with the previous result by the E866 collaboration |4| 
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Table 1 

HERA-B preliminary results on the average transverse momentum (pj) for J/y — ► e + e~ and J/ty — > decays, 
compared to results of previous experiments. 



Target 


Experiment 


Max. pj 


(p T )[GeV/ C ](e+e-) 


( PT )[GeV/ C ]( Al + Al -) 


C, 920 GeV 
W, 920 GeV 
Si, 800 GeV 
Au, 800 GeV 


HERA-B Preliminary 
HERA-B Preliminary 

E771 O 

E789 


4.8GeV/c 
4.8 GeV /C 
3.5GeV/c 
2.6 GeV /C 


1.22±0.01(stat.) 
1.29±0.01(stat.) 


1.22±0.01(stat.) 
l,30±0.01(stat.) 
1.20±0.01 
1.290 ±0.009 
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Figure 4. Nuclear suppression parameter a(xp) for 
two different combinations of carbon and tungsten 
wires, compared to the result of E866 [4 1. 



in the overlap region and allow the measurement to be 
extended to xp = —0.3. 

4.4. y/(2S) Production 

The production cross section for y(2S) mesons is 
measured relative to the J/ty cross section in order to 
reduce systematic uncertainties. In order to extract 



R = 



BR(y(2S) -> £+£-) ■ a v{2s) _ N ¥{2S) £ J/y 

£ V(2S) 



(4) 



the relative yield Ny(2s)/Nj/ V is corrected by the effi- 
ciency ratio Ey v / £y(2s)> which is obtained from a MC 
simulation. 

A preliminary result for the e + e~ channel is R = 
0.13 ±0.02(stat.). As shown in Fig. |5] this result 
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Figure 5. y/"(2S)-to-J/ty production ratio as a function 
of the center-of-mass energy y/s (top) and the atomic 
mass number A (bottom). The preliminary HERA-B 
result (square) is compared to results of previous ex- 
periments. 

agrees well with previous measurements, both as a 
function of the center-of-mass energy and the atomic 
mass number. 

4.5. Xc Production 

The fraction of J/ty produced via the radiative de- 
cay Xcl,2 -> J/Vr -> M + M~7, 



R(Xo) : 



Xc 



(5) 



provides a good test of charmonium production mod- 
els. As depicted in Fig. [6] the Xc signal is extracted 
from the mass difference distribution of the /i + /i~y 
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Figure 6. Xc signal in mass difference distribution 
M([i + [i~y) — M(fi + fi~). Combinatorial background 
has been subtracted. 



and the [i + fi~ systems. The efficiency e r to detect the 
additional photon is determined from a MC simula- 
tion. 

A preliminary analysis of a part of the 2002/2003 
data yields R(x c ) = 0.21 ±0.5(stat.). This result is 
consistent with the previous HERA-B result, R(Xc) = 
0.32±0.6(stat.)±0.4(syst.) 0. 



5. Open Charm Production 

5.1. D Meson Production Cross Sections 

D meson production is studied using minimum bias 
triggered data. The signals shown in Fig. Q are ob- 
tained imposing cuts on the separation of the decay 
vertex from the primary vertex and the impact param- 
eter to the target wire. 

Using these data, the production cross section for 
D° and D + , <r D 0| D + = (a D o + a D +)/2, is determined, 
and the production ratios a D + / cr D o and <r D *+ / <r D o are 
extracted. The cross sections obtained in the lim- 
ited acceptance of the HERA-B detector are extrap- 
olated to the full xp range using Pythia. Prelimi- 
nary results of these measurements are summarized 
in Table |3 As shown in Fig. [8] the result on <r D 0| D + 
agrees well with other measurements at similar beam 
energies. The HERA-B measurement of the ratio 
°b+ / °b° i s m good agreement with pion-induced re- 
sults at lower center-of-mass energies. HERA-B will 
improve the results of previous proton-induced exper- 
iments at similar energies. 
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Figure 8. D° and D + production cross sections (top) 
and cross section ratio (bottom), compared to results 
of previous experiments. Pion-induced results are de- 
picted as open circles, while closed circles stand for 
proton-induced results. The line indicates a prediction 
for the production ratio based on an isospin model 1 6 1 . 



5.2. Limit on D° ^ 

An upper limit on the branching fraction of the 
flavor-changing neutral current process D° — > /i + /x _ 
has been determined analyzing data collected with 
the dimuon trigger. As shown in Fig.|9j three events 
are observed in the signal region, while the expected 
background amounts to 6.0 ± 1.2 events. From these 
numbers, an upper limit for the D° — > branch- 
ing fraction of 2.0 x 10~ 6 is obtained at 90% confi- 
dence level |7|. 

6. Summary 

Measuring the properties of open and hidden charm 
production in proton-nucleus interactions is one of 
the main physics goals of the HERA-B experiment. 
This paper presents preliminary results on a variety 
of physics studies, including open charm and charmo- 
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Figure 7. Open charm decays observed in the minimum bias data-set: D° — > K ;r + (left), D + — > K + w + w (middle), 
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and D* 



D tv — > K 7T+7T+ (right). The corresponding charge-conjugated channels are included. 



Table 2 

Production cross sections and cross section ratios for 
D°, D + , and D* + decays. The cross sections are given 
with statistical and systematic uncertainties. 



Preliminary 


-0.1<x F <0.05 


Full xp range 


<j d o [,ub/nucl] 


21.4±3.2±3.6 


56.3±8.5±9.5 


ct d + [^b/nucl] 


11.5±1.7±2.2 


30.2±4.5±5.8 


Ob«+ [/ib/nucl] 


10.0±1.9±1.4 


27.8±5.2±3.9 


C7 D + / CT D o 




0.54±0.11±0.14 






0.54±0.12±0.10 



nium production cross sections in proton-nucleus col- 
lusions, differential distributions and nuclear effects. 
The HERA-B collaboration is looking forward to fi- 
nalizing these promising analysis efforts in the near 
future. 
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